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A generalized contour deformation method (GCDM) which combines complex rotation and trans- 
lation in momentum space, is discussed. GCDM gives accurate results for bound, virtual (anti- 
bound), resonant and scattering states starting with a realistic nucleon-nucleon interaction. It pro- 
vides a basis for full off-shell t-matrix calculations both for real and complex input energies. GCDM 
competes favorably with analytic continuation. Results for both spectral structures and scattering 
amplitudes compare perfectly well with exact values for the separable Yamaguchi potential. Ac- 
curate calculation of virtual states in the Malfliet-Tjon and the realistic CD-Bonn nucleon-nucleon 
interactions are presented. 

GCDM is also a promising method for the computation of in-medium properties such as the 
resummation of particle-particle and particle-hole diagrams in infinite nuclear matter. Implications 
for in-medium scattering are discussed. 



> 

cn 

o 
m 
o 
m 
o 



X 



I. INTRODUCTION 

The study of two-body resonant structures has a long 
history in theoretical physics, and there exists a variety 
of methods, described in textbooks such as Among 
the more popular methods we have the complex scaling 
method (CSM) and the method based on analytic con- 
tinuation in the coupling constant (ACCC). 

In this work we consider a new approach formulated 
for integral equations in momentum space. The method 
is based on deforming the contour integrals in momen- 
tum space, known as the contour deformation or distor- 
tion method (CDM). It has been shown in Ref. Q that 
a contour rotation in momentum space is equivalent to 
a rotation of the corresponding differential equation in 
coordinate space. The coordinate space analog is often 
referred to as the dilation group transformation, or com- 
plex scaling. The dilation group transformation was first 
discussed and formulated in Refs. |^|^, and was devel- 
oped to examine the spectrum of the Green's function on 
the second energy sheet. 

Complex scaling in coordinate space has for a long 
time been used extensively in atomic and molecular 
physics, see Ref. 0. During the last decade it has also 
been applied in nuclear physics, as interest in loosely 
bound nuclear halo systems has grown, see for example 
Refs. Hj I3- Complex scaling in coordinate space is 
usually based on a variational method and an opti- 
mal variational basis and scaling parameters have to be 
searched for. One of the disadvantages of the coordinate 
space approach is that the boundary conditions have to 
be built into the equations, and convergence may be slow 
if the basis does not mirror the physical outgoing bound- 
ary conditions well. 



There are several advantages in considering the con- 
tour deformation method in momentum space. First, 
most realistic potentials derived from field theoretical 
considerations are given explicitly in momentum space. 
Secondly, the boundary conditions are automatically 
built into the integral equations. Moreover, the Gamow 
states in momentum space are non-oscillating and 
rapidly decreasing, even for Gamow states with large 
widths far from the real energy axis, as opposed to 
the complex scaled coordinate space counterpart. These 
states are represented by strongly oscillating and expo- 
nentially decaying functions. Finally, numerical proce- 
dures are often easier to implement and check. Conver- 
gence is easily obtained by just increasing the number of 
integration points in the numerical integration. 

The contour deformation method (CDM) formulated 
in momentum space is not new in nuclear physics. It 
was studied and a ppl ied in the 1960's and 1970's, see 
for example Refs. [wL Ull , especially in the field of 
three-body systems. These references applied a contour 
rotation method in momentum space. By restricting one- 
self to a rotated contour certain limitations and restric- 
tions however appear in the equations, determined by 
the analytical structure of the integral kernels and poten- 
tials. We will study an alternative approach, by consid- 
ering an extended deformation of the integration contour 
based on rotation followed by translation in the com- 
plex momentum plane. This choice of contour can be re- 
garded as a special case of the Berggren class of contours 
|T3|. Berggren and later Lind studied various 
completeness relations derived by analytic continuation 
of the Newton completeness relation jlj] to the complex 
plane. The Berggren completeness includes discrete sum- 
mation over resonant as well as bound states. Our choice 
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of contour differ from the Berggren class of contours in 
that the contour approaches infinity along complex rays 
in the complex fc-plane as opposed to the various con- 
tours studied by Lind |l4j which approach infinity along 
the real fc-axis. By transforming the momentum space 
Schrodinger equation onto the rotated followed by trans- 
lated contour, we will show that we are able to expose 
and explore more of the physical interesting area on the 
second energy sheet, i.e., the choice of contour enables us 
to study both bound, virtual and resonant states. If one 
restricts the deformation to a rotation of the contour, as 
studied in Refs. [H [H [Hill [3 , we are not able to 
expose virtual states in the energy spectrum, since the 
maximum rotation angle does not allow rotation into the 
third quadrant of the complex momentum plane. This 
limitation is sometimes used as an argument for advo- 
cating different approaches, such as the ACCC method, 
see the recent work of Aoyama jrj. By distorting the 
contour by rotation and translation we are able to in- 
troduce a new feature to the complex scaling method, 
namely accurate calculation of virtual states as well as 
bound and resonant states. Our method represents also 
an alternative to the exterior complex scaling method. 
The exterior complex scaling method was formulated to 
avoid intrinsic non-analyticities of the potential, and in 
this way calculation of resonances in non-dilation ana- 
lytic potentials are possible, see Ref. ^ and references 
therein. By the rotated and translated contour choice 
the CDM is a preferable method compared to all other 
known methods, such as the ACCC method. 

The contour deformation method has also been ap- 
plied to the solution of the full off-shell scattering am- 
plitude (t-matrix), see Refs. 0, El 12, 18]. By rotat- 
ing the integration contour, an integral equation was ob- 
tained with a compact integral kernel. This has numer- 
ical advantages as the kernel is no longer singular. As 
discussed in Ref. [l3|, a rotation of the contour gives 
certain restrictions on the rotation angle and maximum 
incoming/outgoing momentum in the scattering ampli- 
tude. We will again show that our extended choice of 
contour in momentum space avoids all these limitations 
and that an accurate calculation of the scattering ampli- 
tude can be obtained. 

Thus, the method we will advocate allows us to give an 
accurate calculation of the full energy spectrum. More- 
over, it yields a powerful method for calculating the full 
off-shell complex scattering amplitude (i- matrix). It is 
also rather straightforward to extend this scheme to in- 
medium scattering in e.g., infinite nuclear matter. 

In section ^] we outline the contour deformation 
method in momentum space, and in section ITTll the spec- 
tral representation of the full off-shell t-matrix is given 
along with the deformed integration contour. Section 
IIVI presents as a test case a simple separable interac- 
tion which admits analytical solutions for both the en- 
ergy spectrum and the t- matrix. Section Ivl gives numer- 
ical calculations of energy spectra and the t-matrix by 
the contour deformation method. We present numeri- 



cal results for the Yamaguchi T^, Malfliet-Tjon "2^ and 
the charge-dependent Bonn (CD-Bonn) 21] interactions. 
Calculation of virtual states of the CD-Bonn interaction 
is not known by us to have been performed previously. 

II. THEORETICAL FRAMEWORK 

We will in the following use natural units h = c = 1. 
The two-body momentum space Schrodinger equation in 
a partial wave decomposition reads 

TT'^riiik) + - 1 dqq^Vi{k,q)ipni{q) ^ Enii>ni{k). (1) 

2^ TT Jo 

For the sake of simplicity we assume here that the inter- 
action is spherically symmetric and local in coordinate 
space and without tensor components and/or spin-orbit 
coupling. When solving the correspondning equations for 
a more realistic nucleon-nucleon interaction below, these 
degrees of freedom will be accounted for. 

The Fourier-Bessel transform of the potential Vi (r) in 
coordinate space is given by 

/>oo 

Vi{k,k') = / drr^ji{kr)ji{k'r)Vi{r). (2) 
Jo 

The momentum space Schrodinger equation in Eq. 
(with real momenta) corresponds to a hermitian Hamil- 
tonian. The eigenvalues will in this case always be real, 
corresponding to discrete bound states {Eni < 0) and 
a continuum of scattering states {Eni > 0). The eigen- 
states form a complete set, and for a given partial wave I 
the completeness relation, moreprecisely known as reso- 
lution of unity, can be written |l| 

l = El^"')(^»'l + 9 / dk\Mk)){Mk)l (3) 

The infinite space spanned by this basis is given by 
all square integrable functions on the real energy axis, 
known as the space, which forms a Hilbert space. 
Resonant and virtual states can never be obtained by 
directly solving Eq. (Q, as it stands. In a sense, one can 
say that the spectrum of a hermitian Hamiltonian does 
not display all information about the physical system. 

In this paper we study and explore the resonant and 
virtual state spectra by the contour deformation method. 
This is essentially a transformation of Eq. into the 
lower-half complex /c-plane. Such a transformation of 
Eq. can be obtained by an analytic continuation of 
the completeness relation of Eq. ||2Jl to the complex fc- 
plane. One can consider the integral in Eq. Q as an 
integral over the contour T = S -\- C, where the contour 
C is defined on the real fc-axis from — oo to -\-oo and the 
contour S is given by an infinite semicircle in the upper- 
half complex fc-plane closing the contour F. In Ref. 0] 
completeness relations for various inversion symmetric 
contours in the complex fc-plane were derived and dis- 
cussed. Inversion symmetric contours are defined by the 
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following: if z is on C, then —z is also on C. The deriva- 
tion given in Ref. |l4j was based on analytic continuation 
by deforming the contour C defined on the real fc-axis. 
These completeness relations can be regarded as a g ener- 
alization of the Berggren completeness relation [iJl- We 
will therefore label the inversion symmetric contours dis- 
cussed in the extended Berggren class of contours. 
The completeness relation of Berggren was an extension 
of the completeness relation of Eq. Q through the inclu- 
sion of a finite sum over discrete resonant states. By re- 
defining the completeness relation on distorted contours 
in the complex k plane, one can show by using Cauchy's 
residue theorem that the summation over discrete states 
will in general include bound, virtual and resonant states 
|14| . The eigenf unctions will form a biorthogonal set, 
and the normalization follows the generalized c-product 

{{i^nlltpn'l}) = i^nMn'l) = Sn,n' ■ (4) 

The most general completeness relation on an arbitrary 
inversion symmetric contour C = C+ + can then be 
written as 

1 = E / dz\4>i){i;tl (5) 

where C"*" is the distortion of the positive real fc-axis, 
and C~ is the distortion of the negative real fc-axis. The 
symmetry of the integrand has been taken into account, 
that is 

/ dz\i;i){ri\ = I dz\i^i){ri\. 

JC- JC+ 

The summation is over all discrete states (bound, vir- 
tual and resonant states) located in the domain C, de- 
fined as the area above the contour C, and the integral 
is over the non-resonant complex continuum defined on 
C+. The space spanned by the basis given in Eq. (0 
includes all square integrable functions defined in the 
domain C. The complete basis could then be used to 
expand resonant and virtual states defined in the re- 
gion above the distorted contour. Such a complete basis 
is more flexible than a complete basis defined for only 
real energies. From the general completeness relation 
one can deduce the corresponding eigenvalue prob- 
lem, iJlV") — E\ip). The Hamilitonian will in this case 
be complex and non-hermitian, as Gamow and virtual 
states are included in the spectrum. 

In close analogy with the above discussion on complete- 
ness relations, the momentum space Schrodinger equa- 
tion, see Eq. J^l, defined on the positive real fc-axis can 
be continued to the lower-half complex fc-plane. Eq. 
is a Fredholm integral equation of the second kind. Con- 
tinuing Eq. to the lower-half fc-plane, the general rule 
that the moving singularities of the integral kernel must 
not intercept the integration contour must be obeyed, see 
for example Ref. 2] . The choice of distorted contour will 



for each partial wave be based on the a posteriori knowl- 
edge of poles in the scattering matrix. By considering 
the transformed momentum space version of Eq. , the 
choice of contour will in addition be dictated by the an- 
alytic structure of the integral kernel and potential. The 
contour must be chosen in such a way that singularities 
in the potential are located outside the closed integration 
contour. 

In the following we study the analytic continuation of 
Eq. 1^ onto two distorted contours Cf and C^. These 
contours can be regarded as a special case of the extended 
Berggren class of contours. The contour is obtained 
by a phase transformation (rotation) into the lower-half 
complex fc-plane while the second contour C2 will be 
based on rotation followed by translation in the lower-half 
complex fc-plane. These contours differ in an important 
aspect from the extended Berggren class of contours, in 
that they approach infinity along complex rays, and not 
along the real fc-axis. It has previously been assumed as a 
requirement for the choice of distorted contours that they 
approach infinity along the real fc-axis, see for example 
Ref. 113. 

First we consider a contour given by two line seg- 
ments Li and L2. Line Li is given by zi = fcexp {—iO) 
where fc e [0,kmax], L2 by Z2 = kmax ex.p {-i9) (k stih 
real). One can easily show that for an exponentially 
bounded potential in coordinate space the integral in 
Eq. ^ along the arc L2 will go to zero for kmax ^ 00. In 
this case the contour reduces to the line Li. Fig. 
shows the contour along with the exposed and ex- 
cluded two-body spectrum in the complex fc-plane, which 
this contour choice implies. The discrete spectrum con- 
sisting of bound, virtual and resonant states corresponds 
to poles of the scattering matrix Si{k) in the complex 
fc-plane. The contour Ci is part of the inversion sym- 
metric contour Ci = + , as can be seen in Fig. ^ 
For a potential which is analytic in the region above the 
contour and below the real fc-axis, we can then derive 
the transformed Eq. ^ on contour 

(^exp (-2^61)^ - Eni^ ijni{z) + 
2 

exp{-3ie)- dqq^Vi{z,z')i>ni{z') = {), (6) 
TT J a 

where z = fcexp (—16*) and z' = gexp(— J0). Eq. © 
is the momentum space version of the complex scaled 
Schrodinger equation in coordinate space, discussed in 
e.g., Ref. 3=|. A rotation in momentum space, fcexp {—i0), 
is equivalent to the complex scaling r exp {i9) in coordi- 
nate space. The phase-transformation k ^ k exp {~iO) 
is formally a similarity transformation, see for example 
Ref. ,22i]. The restriction on the rotation angle 9 for the 
phase-transformation k k exp {—iO) is given by the re- 
gion of analyticity of the potential. It has been shown in 
Refs. 0,0 that such a transformation does not alter the 
location of bound states in the system; the bound state 
spectrum is invariant under such transformations. This 
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FIG. 1: Contour = Li + L2 is given by the solid line, 
while the contour Cj~ is given by the dashed line. The contour 
Ci = Ci + Ci is clearly inversion symmetric. The two body 
spectrum which is exposed by this contour is marked by filled 
circles • and the excluded spectrum by open circles o. The 
full spectrum includes bound states (BS), virtual (VS), decay 
(DS) and capture (CS) resonant states. 

is also clear from the discussion above on completeness 
relations. The continuum is shifted into the lower-half 
/c-plane, while resonant states will occur as long as they 



are located above the integration contour. The trans- 
formed Eq. © represents a non-hermitian Hamiltonian. 
The eigenvalues are therefore in general complex, and 
the corresponding eigenfunctions form a biorthogonal set 
with the normalization condition given by Eq. Q). 

Next we consider the contour obtained by rotation fol- 
lowed by translation in the lower-half complex /c-plane. 
The contour C2 consists of three line segments. The line 
segment Li is given by a rotation zi — ki exp {—iO) where 
ki G [0, b], L2 is given by a translation Z2 — k2 — ibsm(6) 
where fc2 G [6cos(6'), kmax] and b determines the transla- 
tion into the lower-half fc-plane and L3 by Z3 = kmax — ic 
where c S [6sin(6'), 0]. For kmax 00 the contribution 
to the integral in Eq. along the line segment L3 will 
vanish, and the contour reduces to the line segments 
Li and L2. Fig. [21 shows the contour = Li + L2 + L3 
along with the exposed and excluded two-body spectrum 
which this contour choice implies. The contour C2 is 
part of the inversion symmetric contour C2 — C2 + C2 
clearly seen in the figure. We can now transform the 
momentum space Schrodinger equation onto the dis- 
torted contour given by the lines Li and L2, i.e., 
k,k' z,z'. The integral in Eq. ^ will in this case 
couple the complex momenta zi and Z2- The coupled 
transformed Schrodinger equation then reads 



_]_fzl \ f Mz,)\ 2( J^^dz\z'lVi{z,,z\) JlJz'2z'IVi{z,,z'2)\ f i;^i{z\)\ ^ fMzi)\ 

2^l\0 Zl ) \ i>^l{z2)y^ y f^^dz\z'lVliz2,z'i) J^^dz'2z'lVliz2,z'2) ) \ ^z' 2) ) ^ V'«i(^2) ) ' 

1 



This equation is again a non-hermitian Hamiltonian. The 
basis of eigenstates forms a biorthogonal set, and the nor- 
malization is again given by Eq. Q). The completeness 
relation given in Eq. © , will include a discrete sum over 
bound, virtual and resonant states, and the integration 
along L = Li + L2 is over the complex non-resonant en- 
ergy continuum. 

Fig. O shows a plot of the energy spectra of the phase- 
transformed Eq. © and the rotated + translated Eq. {Tj) 
in the complex energy plane, for a two-body potential 
supporting bound and resonant states. 

Whether one chooses to solve the Schrodinger equation 
on the contour or on the contour C2 depends on 
the problem under consideration. For potentials which 
are analytic in the entire lower-half fc-plane, solving 
along contour is numerically the most straightfor- 
ward method. In most cases the potential has singulari- 
ties in the lower-half fc-plane, and solving on the contour 
C2 enables us to avoid the singularities of the potential, 
while still being able to study resonant structures in the 
system. Another advantage is that there is no restriction 
on the rotation angle 9 as long as the contour is chosen 
so that the poles of the potential are located outside the 



contour. If a potential is of such an analytic structure 
that we are allowed to rotate into the third quadrant of 
the complex fc-plane (see Fig. virtual states will ap- 
pear in the calculated spectra, as long as the potential 
supports virtual states. By solving the Schrodinger equa- 
tion on the distorted contour rotated into the third 
quadrant of the complex fc-plane, we expose a part of the 
negative imaginary fc-axis where virtual states may be 
located, while at the same time excluding a part of the 
positive imaginary fc-axis where bound states may be lo- 
cated. This reminds us that the contour should be chosen 
relative to the partial wave component under study. This 
means that a separate analysis has to be made for each 
partial wave. 



III. TWO-BODY t-MATRIX 

In this section we will discuss different procedures for 
solving the full off-shell t-matrix, and hence the full two- 
body scattering problem. We will consider the general 
mathematical case where the incoming energy is allowed 
to take non-physical values, i.e., the input energy is com- 
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plex. This has relevance for nuclear medium studies 
where the input energy is in general complex. Two meth- 
ods for solving the full ofF-shell t-matrix for arbitrary 
complex energy lo are outlined. 

The t-matrix is defined in operator form by 



t{uj) = V + VG{uj)V, 



(8) 



or 



FIG. 2: Contour — Li + L2 + L3 is given by the solid line, 
while the contour Cj~ is given by the dashed line. The contour 
Ci = Ci + Ci is clearly inversion symmetric. The two body 
spectrum which is exposed by this contour is marked by filled 
circles • and the excluded spectrum by open circles o. The 
full spectrum includes bound states (BS), virtual (VS), decay 
(DS) and capture (CS) resonant states. 
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FIG. 3; General energy spectra in the complex energy plane 
for the complex transformed Schrodinger Eqs. © and Q, 
Fig.a and Fig.b, respectively, bound states are denoted by 
(BS), decay resonant states by (DS) and the complex contin- 
uum by (CC). 



t{uj) =V + VGa{Lu)t{uj). 



(9) 



Here cu is the incoming energy, G'(w) is the resolvent, 
commonly known as the Green's operator, and Go(w) the 
corresponding free Green's operator. In operator form 
they are defined by 



LO - Ha' 
1 

uj-H' 



(10) 
(11) 



The term is the kinetic energy operator and H the full 
two-body Hamiltonian. By expanding the unit operator 
on a complete set of physical eigenstates of H given in 
Eq. we can write the Green's operator as 



G{oj) = 



LO - Eb 



dE, 



l^c)(V'c 
UJ - Ec 



(12) 



Imk 
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FIG. 4: Plot of integration contour C2 where the potential 
singularities are located in the marked area in Fig. a. The ro- 
tation angle 9 is greater than 7r/2 and the integration contour 
encloses a virtual state (A) in the fc-plane. The corresponding 
energy spectrum in the complex energy plane in Fig. b illus- 
trates how the virtual state (VS) is exposed in the spectrum 
and the complex non- resonant continuum (CC). 



This is the spectral decomposition of the Greens's func- 
tion. Here h denotes the discrete bound state spectrum 
and c the positive energy continuum. By projecting 
t{oj) on momentum states, and decomposing into par- 
tial waves, the f-matrix elements ti{k,k']Lo) can be ex- 
pressed as 1-dimensional integral equations. Depending 
on whether we start from JSjl or (jS)) we get the spectral 
or the Fredholm representation of the t-matrix. The con- 
tour deformation method will be applied to the spectral 
representation of the i-matrix while the solution of the 
Fredholm representation of the i-matrix will be based on 
the standard principal value prescription. 

By using Eq. © and Eq. (|12ll we get the spectral rep- 
resentation of the t-matrix by inserting the expansion of 
the Green's function on a complete set of states (|12|l . 
giving 
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Uik,k';co)^mk,k') + ^Yl r dq' q\''v,{k,q)^I^^^^^^^Vi{q',k'). (13) 



The sum over a implies a discrete sum over bound states 
and an integration over the positive energy continuum. 
Eq. is analytic along the real energy axis, except 
for poles located at bound state energies and a branch 
cut along the positive energy axis. In physical two-body 
scattering the incoming energy is defined on the positive 
real energy axis. In this case the integrand in Eq. Ijl^fl is 
singular, and a numerical solution of Eq. H13() is highly 
non-trivial. However, for negative real input energies it 
can be solved by standard numerical procedures. 

We are however interested in the i-matrix for arbi- 
trary complex energies, since an obvious extension of this 
work is to consider in-medium scattering of two nuclcons 



Here C is either or C2 ■ This applies to C" as 
well. We have expanded the Green's operator (|ll|l on 
the basis given in Eq. (O, this implies that a represents 
discrete bound, virtual and resonant states and a non- 
resonant complex energy continuum. From the mathe- 
matical structure of Eq. (|14|l it is easily seen that the t- 
matrix is analytic in lu in the entire complex energy-plane 
except for poles located at bound, virtual and resonant 
state energies. In addition there could also be singular- 
ities in the potential itself. The numerical method for 
solving Eq. H14() is based on matrix diagonalization. We 
first have to diagonalize the corresponding Schrodinger 
equation, and then expand the Green's function on the 
discretized basis obtained. 

Applying this method enables us to obtain ti{k, k';uj) 
for both real and complex energies lu. The integral be- 
comes non-singular on the deformed contour for real and 
positive input energies lo, resulting in numerically stable 
solutions for physical two-body scattering. Eq. H14|) can 
also be considered as an analytic continuation of Eq. H13|) 
for complex input energies u, and stable numerical solu- 
tion can be obtained for complex energies above the dis- 
torted contour. The limitation of this method is due to 
that most potentials in momentum space have singulari- 
ties in the complex plane when one argument is real and 
the other is complex. By applying contour , which is 
based on rotation into the complex plane, in most cases 
there will be restrictions on both rotation angle {9) and 
maximum incoming and outgoing momentum (k,k'), see 
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or to study the resummation of large classes of many- 
body diagrams. In a nuclear medium calculation the 
self-consistently determined quasiparticle energies are in 
general complex. 

We can achieve this by an analytic continuation of the 
contour integrals into the complex /c-planc. This rep- 
resents a deformation of the integration contour in the 
complex fc-plane. See Refs. for validity and math- 

ematical proofs concerning analytic continuation of inte- 
gral equations. 

We will consider the solution of Eq. (|13|l by integrating 
over the two contours, or C^, discussed in Sec. m 
resulting in the following expression for i/(fc, k'; lu) 



(14) 

I 

for example Ref. [l]| . 

Using contour we can avoid these limitations by 
choosing the integration contour in such a way that the 
potential singularities always will lie outside the integra- 
tion contour, and therefore do not give any restriction 
on rotation angle and maximum incoming and outgoing 
momentum. 

The partial wave decomposition of Eq. Q gives 
the Fredholm representation, commonly known as the 
Lippmann-Schwinger equation 

Jo E{q) 

In physical two-body scattering the input energy is real 
and positive. In this case the Fredholm integral Eq. 115|) 
has a singular kernel of Cauchy type. Solving singular in- 
tegrals can be done by Cauchy's Residue theorem, where 
we integrate over a closed contour enclosing the poles. 
There are two ways of doing this, either by letting z lie an 
infinitesimal distance above the real axis, i.e., lu — *■ z-\-ie, 
or by letting lo lie on the real axis. In both cases we must 
choose a suitable contour enclosing the singularity. If we 
choose the latter position of the singularity, we get a 
Cauchy principal-value integral where we integrate up to 
- but not through - the singularity, and a second contour 
integral, where the contour can be chosen as a semicir- 
cle around the singularity. Eq. H15|) can thus be given in 
terms of a principal value part and a second term com- 
ing from integration over the semicircle around the pole. 
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ti{k,k',Lu) = Viik,k') + ^y f f dzdz' z^z'\{k,z)^^^^^^^^^Vi{z\k% 

„ JC Jc LU- Ea 
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The result is 



ti{k,k';uj) = Vi{k,k') 



-V 



dqq^Vlik,q)ti{q, k') 
LO - E{q) 



2ifikoVi{k, ko)ti{ko, k';uj). 



(16) 



r 



By rewriting the principal value integral using the rela- 
tion 



V 



dk 



fjk) 
Aiq — k^ 



we obtain an equation suitable for numerical evaluation. 
Eq. I|16|l can be converted into a set of linear equations 
by approximating the integral by a sum over N Gaussian 
quadrature points (fc^; j = 1, ...N), each weighted by Wj 
The full off-shell t-matrix is then obtained by matrix in- 
version. This method for solving integral equations is 
known as the Nystrom method. It is numerically effec- 
tive and stable, except for the rare case when the incom- 
ing energy u coincides with or is very close to one of the 
integration points. 

So far we have only considered physical input ener- 
gies in Eq. (|15|l . but it has been shown in Ref. '5| that 
the analytically continued Lippmann-Schwinger equation 
to complex energies takes the same form as Eq. H16(l . By 
solving the full off-shell t-matrix for arbitrary complex in- 
put energy, we do not have to alter the set of linear equa- 
tions obtained for physical energy, the only modification 
is that the energy is complex. It should be noted that the 
Lippmann-Schwinger equation H14|l can be solved by the 
contour deformation method as well, giving a compact in- 
tegral kernel for positive incoming energies. By contour 
distortion the principal value prescription is avoided, and 
a numerical solution of Eq. (|14|l is stable even for incom- 
ing energies coinciding with the integration points. 



IV. AN ANALYTICALLY SOLVABLE 
TWO-BODY POTENTIAL 

We now consider a separable potential given by Yam- 
aguchi pjj . It models s- and p- waves. A separable in- 
teraction in momentum space is analytically solvable, see 
for example Ref. |3] for a demonstration. The Yamaguchi 
interaction therefore admits analytical solution of the full 
off-shell i-matrix and the f-matrix poles, corresponding 
to the energy spectrum. The Yamaguchi s-wave potential 
supports bound and virtual states, while the Yamaguchi 
p-wave potential supports bound, virtual and resonant 
states. The Yamaguchi potential is therefore very use- 
ful in modelling loosely bound two-body systems which 
may have a rich resonant state structure, and for check- 
ing the numerics in our calculations of the i-matrix and 
the energy spectrum. 



The s-wave Yamaguchi potential has the form 

Fo(fc,'z) = -A5o(%o(g), (18) 

where 

1 



(19) 



In natural units and mass 2^ = 1 MeV, the full off-shell 
t-matrix for the s-wave potential reads P| 



to{k,q;z) = -A- 



9o{k)ga{q) 



l-XlJ-dk^ 



(20) 



-E (k^+p-^) 

The integral in the denominator can be evaluated, giving 

1 _X{(3 + iVE)^ 



9 r°° 

A- / dk ——— - 

TT Jo k^-E{k^+ (3^) 2 /3(/32 + Ef 

The s-wave i-matrix in closed form reads 

.9o(fc)go(g) 



h{k,q-E) - -A- 



1 - - 



(21) 



(22) 



Writing k = —ik where E = k^ we can solve for the 
t-matrix poles as zeroes of the denominator 



1 



Solving for k gives 



A (/3 - k)^ 
2 /3(/32 _ k2)2 



-f3± 



0. 




(23) 



(24) 



and we see that the poles are located along the imagi- 
nary fc-axis. Poles of the ^-matrix located on the positive 
imaginary k-eods represents bound states, while poles lo- 
cated on the negative imaginary fc-axis represents virtual 
states, giving a bound (k > 0) and a virtual state (k < 0) 
for A > 2/3^ and for A < 2f]^ two virtual states. 
The separable p-wave interaction is given by 



where 



Viik,q) = -Agi(fc)gi(g), 



ffi(fc) 



fc2 + ^2 ' 



(25) 



(26) 
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and the t-matrix now becomes 0] 

9i{k)gi{q) 



ti{k,q-E) = -A 



1 >2 f°° ^i, fc!__fci„ 
^ -^7 Jo "■'^ 



(27) 



Solving again the integral in the denominator gives the 
t-matrix in closed form 



ti{k,q-E) = -A 



9i{k)gi{q) 



^ - 2 



and solving for the poles gives in terms of k = —ik 
« = ^-/3±iv/4A2-8A/3. 



(28) 



(29) 



We see that for p-waves the interaction supports bound, 
virtual and resonant states. The bound state condition 
is 



A > 2(3. 



(30) 



giving in addition a virtual state. The interaction has 
a branchpoint at fc = where the bound and virtual 
state meet and move symmetrically from the imaginary 
axis into the lower-half fc-plane giving capture and decay 
resonant states. Fig. |S1 shows the pole trajectory for the 
Yamaguchi p-wave interaction. 



A Im k 



Re k 




FIG. 5: Pole trajectory for the p-wave Yamaguchi interaction 
in the complex fc-plane. 



A. Results for the Yamaguchi potential 

The Yamaguchi potential is analytic in the entire com- 
plex /c-plane except for singularities located at fc, /c' = 
±Z/9, arising from the separable part of the potential 
l/{k'^+0^). This singularity gives a restriction on the in- 
tegration contour; if we consider the rotated contour Cj^, 
the rotation angle 9 has to be less than 7r/2, i.e., 9 < tt/2, 
in the complex fc-plane. This contour choice will not be 
able to give the virtual states in the calculated energy 
spectrum. The contour could on the other hand be 
chosen to lie above the singularity, i.e., 5sin(6') < (3. By 
this choice there is no restriction on the rotation angle 9 
and we can rotate into the third quadrant of the complex 
fc-plane. By this procedure virtual states can be included 
in the calculated energy spectrum. 

In Tabled we present numerical versus exact values for 
the s-wave virtual states in the Yamaguchi potential, us- 
ing contour C^. The contour was rotated 9 = 2Tr/3 
into the third quadrant and translated c = 5sin(27r/3) w 
4.3301 MeV in the lower-half complex fc-plane. Fig. El 
shows the contour along with the excluded spectrum and 
the singularities of the s-wave Yamaguchi potential giv- 
ing restrictions on the contour choice. The parameter 
/3 was held fixed, and equal to (3 = 6 MeV. The poten- 
tial depth A was chosen to give only virtual states. The 
contour choice exposes one virtual state while one virtual 
state will be excluded from the calculated spectrum. The 
calculations used A^l = 50 integration points along line 
Li and N2 = 50 integration points along L2 
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V. NUMERICAL RESULTS 

Below we present numerical calculations for the com- 
plete ei ierg y spectrum of the separable Yamaguchi po- 
tential for the Malfliet-Tjon 20| and the realistic 
charge dependent Bonn (CD-Bonn) 21] nucleon-nucleon 
interactions. We also present calculations of the t-matrix 
for the given interactions within the formalism outlined 
above. The analytic structure of the various interactions 
is of importance for the choice of the integration contour 
in the complex plane, as will be discussed. 



FIG. 6: Excluded (o) and included (•) spectrum for the s- 
wave Yamaguchi potential, with /9 = 6 MeV and A chosen 
to give only virtual states (VS) in the full spectrum. The 
calculation of the spectrum for this contour choice is given in 
Table Potential singularites are given by (x) in the figure. 

In Table ^and Table Ulll we present numerical versus 
exact values of the p-wavc energy spectra, calculated us- 
ing contour and . The calculations used VI = 50 
integration points along line Li and V2 = 50 integra- 
tion points along L2 for contour C^. For the contour 
a number of V = 100 integration points was used. 
The rotation angle is the same for both contours and 
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equal to 9 = 7r/6. The parameter j3 in the potential was 
held fixed, and equal io (5 — & MeV. The translation pa- 
rameter of contour C2 is given hy b = 3.5 MeV. Fig. [3 
shows the contour along with the excluded and included 
spectrum for the contour and the contour C2 ^ along 
with the singularities of the p-wave Yamaguchi potential 
giving restrictions on the contour choice. 

i Imk 

x+P =6 MeV 
O BS 



b=3.5MeV Re k 




=-6MeV 



FIG. 7: Excluded (o) and included (•) spectrum for the p- 
wave Yamaguchi potential, with 13 = 6MeV and A chosen to 
give bound (BS), virtual (VS), capture (CS) and decay (DS) 
states in the full spectrum. The calculations of the spectrum 
for contour {CI) and contour {C2) are given in Ta- 
bles [H] and ^HJ respectively. Potential singularites are given 
by (x) in the figure. 





Exact 


Contour C2 


A 


Real 


Real 


200 


-3.67687178 


-3.67687178 


250 


-2.06107759 


-2.06107759 


300 


-1.00000000 


-1.00000000 


350 


-0.35925969 


-0.35925969 


390 


-0.08947449 


-0.08947449 



TABLE I: Numerical calculations for the Yamaguchi s-wave 
virtual states versus exact values. Energies are given in units 
of MeV. The numerical calculations were obtained using the 
contour deformation method (CDM) along contour C^. 

The results for the virtual states for the p-wave Yam- 
aguchi potential using contour C^, are given in Table Hvl 

We conclude that there is no difference between the nu- 
merically calculated and the exact bound, virtual and res- 
onant energy spectra. To see a difference for the number 
of integration points used in the calculations we have to 
include up two 12 significant digits. Using double preci- 
sion we get for the p-wave (A = 13) virtual state the exact 
value E = -1.697224362268005 MeV while the numeri- 
cal resuh in double precision is E = -1.697224362276472 
MeV. The calculated results for resonant states with large 
widths, i.e., far from the real energy axis agree per- 
fectly with exact values. This illustrates the advantage 



Exact Contour C2 



A 


Real 


Imag. 


Real 


Imag. 


13 


-5.30277586 





-5.30277586 


0. 


12.6 


-2 80486369 





-2 80486369 


0. 


12.2 


-0.77620500 





-0.77620500 


0. 


11.8 


0.57999998 


-0.15362291 


0.57999998 


-0.15362291 


11.7 


0.85500001 


-0.28102490 


0.85500001 


-0.28102490 


11.5 


1.37500000 


-0.59947896 


1.37500000 


-0.59947896 


11.0 


2.50000000 


-1.65831244 


2.50000000 


-1.65831244 



TABLE II: Numerical calculations for the Yamaguchi p-wave 
bound and resonant energy spectra versus exact values. En- 
ergies are given in units of MeV. The numerical calculations 
were obtained using the contour deformation method (CDM) 
along contour C^. 



Exact Contour C 



A 


Real 


Imag. 


Real 


Imag. 


13 


-5.30277586 





-5.30277563 


0. 


12.6 


-2.80486369 





-2.80486362 


0. 


12.2 


-0.77620500 





-0.77620499 


0. 


11.8 


0.57999998 


-0.15362291 


0.57999998 


-0.15362291 


11.7 


0.85500001 


-0.28102490 


0.85500001 


-0.28102490 


11.5 


1.37500000 


-0.59947896 


1.37500000 


-0.59947896 


11.0 


2.50000000 


-1.65831244 


2.50000000 


-1.65831244 



TABLE III: Numerical calculations for the Yamaguchi p-wave 
bound and resonant state energy spectra versus exact values. 
Energies are given in units of MeV. The numerical calcula- 
tions were obtained using the contour deformation method 
(CDM) along contour ■ 



of the contour deformation method in momentum space 
as opposed to the complex scaling analog in coordinate 
space, where convergence of these states will expectedly 
be slower due to the slowly decreasing oscillating wave- 
functions. Fig. |S1 gives a plot of the calculated energy 
spectrum for the p-wave Yamaguchi potential, using con- 
tours Ci and C2 ■ Contour is rotated by 6* = 27r/5 
in the energy plane while contour C2 is rotated into the 
third quadrant of the complex /c-plane hy 6 — 37r/5 giving 





Exact 


Contour 


A 


Real 


Real 


13 


-1.69722438 


-1.69722438 


12.6 


-1.15513635 


-1.15513635 


12.2 


-0.46379500 


-0.46379500 


12.1 


-0.25000000 


-0.25000000 



TABLE IV: Numerical calculations for the Yamaguchi p-wave 
virtual state energy spectra versus exact values. Energies are 
given in units of MeV. The numerical calculations were ob- 
tained using the contour deformation method (CDM) along 
contour with 9 > 7v/2. See Fig.|S]for a plot of the contour 



10 



a rotation oi 9 = 6tt/5 in the energy plane. The resonant 
state appears at the same location for both contours. 



Complex energy plane 



C2 




RS 



2 
1.5 
1 

0.5 


-0.5 

-1 
-1.5 

-2 
-2.5 

-3 



FIG. 8: Calculated energy spectra for the p-wave Yamaguchi 
potential with parameters X — 11, P = 6 using the contours 
(C1),C^ {C2). The discretized non-resonant continuum 
lies along the contours, and a resonant state (RS) appears in 
the energy spectra. 



Next we present a calculation of the t-matrix elements 
ti(k,k;E) for both real and complex input energy (see 
Table |V]and Table lVll ') for the p-wave Yamaguchi poten- 
tial. The ^-matrix is calculated by applying the spectral 
representation on the deformed contour . The contour 
is rotated into the third quadrant of the complex fc-plane 
and the difference between numerical calculations and 
exact values is of the same order as for the calculations 
of the energy spectrum. The rotation, translation and 
number of integration points were the same as above. 





Exact 


Contour 


k 


Real 


Imag. 


Real 


Imag. 


1.5 


0.102303714 


-0.061759732 


0.102303714 


-0.061759732 


3 


0.295657724 


-0.178485632 


0.295657724 


-0.178485632 


4 


0.425747126 


-0.257019311 


0.425747126 


-0.257019311 


6 


0.461965203 


-0.278883785 


0.461965203 


-0.278883785 


7.5 


0.439705133 


-0.265445620 


0.439705133 


-0.265445620 


9 


0.393627137 


-0.237628803 


0.393627137 


-0.237628803 


10 


0.342892975 


-0.207001090 


0.342892975 


-0.207001090 



TABLE V: Calculation of the p- wave {I = 1) i;(fc, k, E) for real 
momenta k and input energy E = 5 MeV. ti{k, k, E) is given 
in units of MeV~'^. The interaction parameters are A = 12.5, 
/3 = 6 MeV, the rotation angle is 6^ = 27r/3 while the trans- 
lation is given by 6 = 3.5 MeV. The number of integration 
points was A^l = Af2 = 50 





Exact 


Contour 


k 


Real 


Imag. 


Real 


Imag. 


1.5 


0.094043918 


-0.030566057 


0.094043918 


-0.030566057 


3 


0.271786928 


-0.088335901 


0.271786928 


-0.088335901 


4.5 


0.391373158 


-0.127203703 


0.391373158 


-0.127203703 


6 


0.424667060 


-0.138024852 


0.424667060 


-0.138024852 


7.5 


0.404204220 


-0.131374046 


0.404204220 


-0.131374046 


9 


0.361846477 


-0.117606975 


0.361846477 


-0.117606975 


10.5 


0.315208495 


-0.102448739 


0.315208495 


-0.102448739 



TABLE VI: Calculation of the p-wave ti{k,k,E) for real 
momenta k and complex input energy E = 5 — 2i MeV. 
ti{k,k,E) is given in units of MeV"'^. Interaction parame- 
ters are A — 12.5, /3 = 6 MeV , the rotation angle is = 27r/3 
while the translation is given by 6 = 3.5 MeV. The number 
of integration points was 7V1 — N2 = 50. 



B. Results for the Malfliet-Tjon interaction 

The Malfliet-Tjon interaction j23| is a superposition of 
Yukawa terms. This interaction resembles the form of a 
realistic nucleon-nucleon interaction with attractive and 
repulsive parts in the ^5*0 channel. It can be fitted to re- 
produce the ^S'o phase shift in nucleon- nucleon scattering 
rather well. The interaction in coordinate representation 
is given by 



V{r) = Va 



exp {-fJ-Ar) 



Vb 



exp (-/isr) 



(31) 



This interaction will support bound and virtual states for 
s-waves, while for higher angular momentum resonances 
will appear. It is known that the channel in the 
nucleon-nucleon interaction supports a virtual state near 
the scattering threshold. 

We will only consider s-waves and calculate the en- 
ergy spectrum using contour rotated into the third 
quadrant of the complex /c-plane. First we consider the 
analytic structure of the Malfliet-Tjon interaction in mo- 
mentum space. The Fourier-Bessel transform of H31|) . for 
the s-wave, is 



Vb 



4kk' 
1 

ifcfc' 



{k + k'f 



/4 



In 



(fc- 
(fc- 



-k'y 



(32) 



Sk-k')^ + ^il^ 

By analytic continuation of the interaction H32|) to the 
complex fc-plane, there will be singularities in the inter- 
action for 



^')'+m1b = 0, 



and 



2 

t^A,B 



0. 



(33) 



(34) 



Eq. H33|) is only satisfied if one of the arguments is real 
and the other complex. For z' = k' real, we have singu- 
larities for ^[z] = =Ffc' A Q[z] = ±fi. Eq. is satisfied 
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for (z + z') = ±ifiA,B, which gives 3?[z] = ^[z'] A = 
^fJ'A,B — '^i^']- Solving the eigenvalue problem by the 
contour deformation method on the purely rotated con- 
tour , singularities appear along the imaginary axis 
given by 5[z] — ±^a,b — '^i^']- This gives a restriction 
on rotation angle, 9 < tt/2. This docs not allow a ro- 
tation into the third quadrant of the complex fc-plane, 
and calculation of virtual states is not possible by this 
contour choice. This problem is resolved by solving the 
eigenvalue problem on the rotated -I- translated contour 
C2 ■ For a rotation angle 9 > Tr/2 singularities appear on 
the contour for 3?[z] = 5R[z'] A + ^[z'] = ±ha,b- 
By imposing a lower bound on the translated line seg- 
ment L2, given by c < ^a,b/2, we avoid all singularities. 
By this choice, rotation into the third quadrant of the 
complex fc-plane is allowed, and an exploration of virtual 
states is possible. 

If we want to extract the t-matrix along the real fc-axis 
by the contour deformation method, see Sec. IIIII singu- 
larities appear when Eq. H33|) is satisfied. If we integrate 
along contour there will always be a singularity on 
the contour given by 



Imk 



where 



in, 



ji/ tan(6'). 



(35) 



(36) 



and /X = minlfiA, fJ's]- For k, k' > kmax the contour 
will pass through the singularity of the interaction and 
Cauchy's integral theorem cannot be applied. However, if 
the interaction, Vi(fc, /c'), is approximately zero for fc, k' > 
kmax integrating along contour can be done as long 
as 



< arctan (/i/fc„ 



(37) 



In this sense we may call kmax the cutoff momentum. 
This choice may cause numerical unstable solutions for 
small values of momenta, since the rotated contour may 
lie very close to the real fc-axis where the integral kernel is 
singular. The same conclusion has already been pointed 
out in Ref. [l^. 

Here we see the advantage of integrating along contour 
C2', Not only will it be able to reproduce virtual states in 
the spectra, but it will also give accurate calculations of 
the i-matrix for real incoming and outgoing momentum. 
It will always be possible to choose a contour C2 lying 
above the nearest singularity z — ^[k] — z/i, implying 
no restriction on rotation angle 9 irrespective of cutoff 
momentum kmax- Fig. (Ogives an illustration. 

A calculation of resonant and virtual states in the 
Malfliet-Tjon interaction is given in Ref. |2^. The t- 
matrix poles were calculated by analytic continuation of 
the t-matrix to the second energy sheet, and thereby 
searching for poles. In Table IVIII a calculation of the 
virtual neutron-proton s-wave states with reduced mass 
M„p = 938.926 MeV and c = 1 for interac- 

tion parameters Va = 7.291 MeV, fiA = 613.69 MeV, 
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FIG. 9: Illustration of potential singularities for the Malfliet- 
Tjon interaction Vi{k,z'), where k is real and z' complex. 
Restriction on the integration contours {CI) and {C2) 
is clearly seen for given cutoff' momentum kmax {Vi{k, k') « 

for k,k > kmax)- 



Hb = 305.86 MeV held fixed while Vb was varied as 
given in Tabic IVTll The rotation is 9 — 27r/3 and the 
translation c = 100 sin(27r/3) MeV. The convergence is il- 
lustrated by increasing the number of integration points. 
We see a small difference in the calculated values for the 







CDM 




Rof.f24] 




A 


B 


C 


D 


Vb 


Energy 


Energy 


Energy 


Energy 


-2.6047 


-0.066674 


-0.066653 


-0.066653 


-0.06663 


-2.5 


-0.310114 


-0.310115 


-0.310115 


-0.31004 


-2.3 


-1.229845 


-1.229845 


-1.229845 


-1.22970 


-2.1 


-2.679069 


-2.678979 


-2.678979 


-2.67878 



TABLE VII: Calculation of the neutron-proton virtual state 
as function of interaction parameter A for the s-wave Malfliet- 
Tjon interaction using the deformed integration contour C^. 
The convergence is illustrated by increasing the number of 
integration points. Column A used Nl = 30, N2 = 50 inte- 
gration points, column B used Nl = 100, N2 = 100 integration 
points and column C used Nl = 150, N2 = 250 integration 
points. Comparison is made with the calculations of Ref. |2^ . 
given in column D. 



virtual state in the Malfliet-Tjon interaction. As we ob- 
tained convergence of the virtual state energies by in- 
creasing the number of integration points, this suggests 
that our results are stable. Fig.llOlshows a plot of the cal- 
culated t-matrix elements ti(k, fc; E) for incoming energy 
E — 100 MeV for the s-wave Malfliet-Tjon interaction, 
using the contours and C^. In this case we used a 
rotation angle 9 = tt/Q for both contours, and a trans- 
lation b = 300 MeV for contour . The Malfliet-Tjon 
interaction has a singularity along contour Cj^, given by 
kmax = m/ tan(6') = 529.77 MeV. In this case the inter- 
action is not sufficiently neglible for k^k' > kmax, and 
using contour will not give accurate calculation of 
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the t-matrix, this is clearly seen in the plot of the cal- 
culated results. This illustrates clearly the advantage of 
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FIG. 10: Plot of f-matrix elements for the s-wave Malfliet- 
Tjon interaction using the contour (dashed line) and the 
contour (solid line). The potential singularity along con- 
tour Cj*" is clearly displayed, and located at kmai. = 529.77 
MeV. 

integrating along the contour instead of . 

A useful check of the numerics is the on-shell unitarity 
for the S'-matrix. Table IVIIII reports calculations done 
by the principal- value prescription and the contour defor- 
mation method C^- We used TV = 50 integration points 
for the principal value integration, for the contour defor- 
mation method we used A^l = 30 and iV2 = 100, 200. 
We observe that we need a higher number of integration 
points along line L2 for obtaining convergence. 



C. Results for the CD-Bonn nucleon-nucleon 
interaction 

Finally, we present a calculation of the energy spec- 
trum of the charge-dependent Bonn interaction (CD- 
Bonn) and phase shifts for a set of selected partial waves. 
The CD-Bonn interaction is given in Ref . |2j . We illus- 
trate the power of the contour deformation method by 
reproducing the deuteron bound state and the virtual 
states in the ^Sq isospin triplet channel. In addition, 
phase shifts are calculated within the same method, and 
the results agree perfectly with other theoretical predic- 
tions. 

The tensor component of the CD-Bonn interaction cou- 
ples angular momentum, Z = j — 1, Z = j -I- 1, in the spin 
triplet channel. It is straightforward to include this cou- 
pling in the formalism outlined in the previous sections. 
Instead of NxN matrices we get 2Nx2N matrices in the 
coupled case, for the uncoupled channels we still have 
NxN matrices. We should of course label the equations 
with quantum numbers L, S, J, T. 



PV CDM 



it [MeV] 




iV = 50 


7V2 = 100 


7V2 = 200 


10. 


1 


00000000 


1.00000811 


1.00000811 


110. 


1 


00000000 


0.999999940 


1.00000000 


210. 


1 


00000000 


0.999999940 


1.00000000 


310. 


1 


00000000 


0.999999881 


1.00000000 


410. 


1 


00000000 


0.999999881 


1.00000000 


510. 


1 


00000000 


0.999999881 


1.00000000 


610. 


1 


00000000 


0.999999881 


1.00000000 


710. 


1 


00000000 


0.999999821 


1.00000000 


810. 


1 


00000000 


0.999999821 


1.00000000 


910. 


1 


00000000 


0.999999821 


1.00000000 


1010. 


1 


00000000 


999999762 


1 00000000 


1110. 


1 


00000000 


0.999999702 


1.00000000 


1210. 


1 


00000000 


0.999999702 


1.00000000 


1310. 


1 


00000000 


0.999999642 


1.00000000 


1410. 


1 


00000000 


0.999999523 


1.00000000 


1510. 


1 


00000000 


0.999999642 


1.00000000 



TABLE VIIL Calculations of 5- matrix norm, 15(^)1, for the 
s-wave Malfliet-Tjon potential with parameters Va ~ 7.291 
MeV, HA = 613.69 MeV, Vb = -5 MeV and /is = 305.86 
MeV. Column 2 gives results for the principal value prescrip- 
tion (PV), while columns 3,4 give results for the contour de- 
formation method (CDM). 

The realistic nucleon-nucleon interaction does not have 
a resonant structure in the low energy region E < 350 
MeV. It does however support virtual states in the ^Sq 
isospin triplet channel and a bound state in the coupled 
isospin singlet channel (the deuteron) '^Si — "^Di. We 
compare the calculated virtual state locations in the com- 
plex fc-plane with the values obtained by the effective 
range approximation, see Ref. ill, fl- In the following 
we do not include Coulomb effects when considering the 
isospin tz — — I channel (proton-proton scattering). 

The effective range approximation for the s-wave poles 
is given by 







1 


rNN ^ 


y rNN\aNN\ 





The theoretical (see below) and experimental values 
for the ^Sq scattering lengths unn and effective range 
rjvAf are given in Table Hxl 

To apply the contour deformation method by inte- 
grating along contour C2 we first have to locate pos- 
sible singularities in the CD-Bonn interaction. The CD- 
Bonn interaction is given explicitly in momentum space. 
The derivation of the interaction is based on field the- 
ory, starting from Lagrangians describing the coupling 
of the various mesons of interest to nucleons The 
one-boson-exchange interaction is proportional to 

y(k,q)cx T4(k,q)F2(k,q;A„). (39) 



13 



CD-Bonn Experiment 



dpp 


-17.4602 






2.845 




dnn 


-18.9680 


-18.9 ±0.4 


fnn 


2.819 


2.75 ±0.11 


Clnp 


-23.7380 


-23.740 ±0.020 


fnp 


2.671 


2.77 ±0.05 



TABLE IX: Scattering lengths (a) and effective ranges (r) 
for the ^5*0 channel, in units of fm. For the proton - proton 
channel Coulomb effects are not included. 

Both T4((k, q) and F^(^, q; Aq.) contain terms of the form 



(k-q)2±m2' ^""^ 

which are of the Yukawa type. These terms determine the 
analyticity region of the interaction in the fc-plane (see 
discussion of the Malfliet-Tjon interaction above. We see 
that the poles of the interaction are determined by the 
various meson masses ttIq and cut-off masses A^. Con- 
sidering the solution of the eigenvalue problem by the 
contour deformation method using contour C^, singu- 
larities in the interaction appear for 9 > 7r/2, i.e., a ro- 
tation into the third quadrant of the complex fc-plane. 
The maximum translation into the complex fc-plane is 
then determined by the smallest meson mass entering 
the potential, which is the 7r-meson, m^o — 134.9764 
MeV. For a given rotation into the third quadrant, i.e., 
9 > 7r/2, we get a restriction on the translation parame- 
ter; b < 134.9764/2 sin( 6*). Here b is given in MeV, as we 
are still using natural units. 

Table IX) gives results for the virtual states in the 
channel by the contour deformation method. A com- 
parison with the effective range calculation of the vir- 
tual state poles is also shown. The contour was rotated 
hy 9 = 27r/3 into the complex fc-plane and translated 
c = -30sin(57r/7) MeV or c « -23.45 MeV in the lower- 
half fc-plane. This is sufficient to reproduce the virtual 
states in the ^5*0 channel, as they are known to lie very 
close to the scattering threshold, k « —hcO.OBi MeV. By 
this contour choice the full energy spectrum is obtain- 
able since it is known a posteriori that the ^Sq channel 
supports only virtual states near the scattering thresh- 
old. The convergence of the numerical calculated val- 
ues is demonstrated by increasing the number of integra- 
tion points N. Table IXll illustrates the convergence of 
the deuteron binding energy as a function of integration 
points. The contour was rotated hy 9 = tt/6 into the 
complex fc-plane and translated c = —/ic sin(7r/6) MeV 
or c « —100 MeV into the lower-half fc-plane. 

Figs. I11I12I and El present calculations of nucleon- 
nucleon phaseshifts for the uncoupled ^Sq isospin triplet 
channel, by the contour deformation method. The con- 
tour was rotated by 27r/3 into the third quadrant of the 
complex fc-plane, and translated c = 50sin(27r/3) MeV. 
We used iVl = 30 and N2 = 70 integration points along 



CDM EFR 





a 


b 


c 






Energy 


Energy 


Energy 


Energy 


-1 (PP) 


-0.11766 


-0.11761 


-0.11761 


-0.11763 


1 (nn) 


-0.10070 


-0.10069 


-0.10069 


-0.10070 


(np) 


-0.06632 


-0.06632 


-0.06632 


-0.06632 



TABLE X: Calculation of virtual state energies in the ^So 
isospin triplet channel by the effective range approximation 
(EFR) and the contour deformation method (CDM), in units 
of MeV. Convergence is obtained in each isospin channel. Col- 
umn (a) used A''l — 20, A'^2 = 30 integration points, column 

(b) used A''l = 20, N2 = 50 integration points while column 

(c) used A''l = 30, A'^2 — 80 integration points 



Nl N2 Energy 
20 30 -2.224581 
20 50 -2.224573 
30 80 -2.224574 
50 100 -2.224575 
50 150 -2.224575 

TABLE XI: Calculation of deuteron binding energy in the 
coupled channel "^Si —^Di by the contour deformation method 
(CDM), in units of MeV. See Ref. ^ for further details. 



line Li and line L2, respectively. Our calculation is given 
by the solid line while the circles are calculations given in 
Ref. [231 . For proton-proton scattering Coulomb effects 
are not included. 
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FIG. 11: Phaseshift for proton-proton scattering in the ^So 
channel. Coulomb effects are not included in the calculations. 
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FIG. 12: Phaseshift for neutron-proton scattering in the ^So 
channel. Calculation by the contour deformation method 
is given by the solid line. Circles are calculations given in 
Ref. m. 
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FIG. 13: Phaseshift for neutron-proton scattering in the ^So 
channel. Calculation by the contour deformation method 
is given by the solid line. Circles are calculations given in 
Ref. 

VI. CONCLUSIONS AND PERSPECTIVES 

A generalized contour deformation method in momen- 
tum space has been presented. The deformation of the 
integration contour is generalized to rotation followed 



by translation in the complex fc-plane. This generaliza- 
tion makes it possible to handle both dilation and non- 
dilation invariant potentials. We have shown this to be 
a powerful procedure for studying resonances and virtual 
states in two-body systems. The method has also been 
successfully applied to the calculation of the full off-shell 
t-matrix. The aim of this work has been to establish the 
formalism for the free scattering case, basing the analysis 
on schematic and realistic nucleon-nucleon interactions. 
Work on extension to few-body Borromean halo systems 
|25l | with more than two constituents, is in progress. 

The exposed formalism allows for stable numerical cal- 
culations of bound states, resonances and virtual states, 
in addition to yielding a fully complex on-shell and off- 
shell scattering matrix, starting with a realistic nucleon- 
nucleon interaction. 

This allows for several interesting applications and ex- 
tensions. As is well-known, one of the major challenges 
in the microscopic description of weakly bound nuclei 
is a proper treatment of both the many-body correla- 
tions and the continuum of positive energies and decay 
channels. Such nuclei pose a tough challenge on tradi- 
tional nuclear structure methods, based essentially on the 
derivation of effective interactions and the nuclear shell- 
model, see for example Ref. for a review. In the 
traditional approaches only bound states typically en- 
ter the determination of an effective interaction, be it 
either based upon various many-body schemes or more 
phenomenologically inspired approaches. Coupled with 
large-scale shell model studies, several properties of sta- 
ble nuclei are well reproduced. However, weakly bound 
nuclei may have a strong coupling to unbound states, 
either resonances or virtual states, as described in for ex- 
ample Refs. [23,12^ 113 • This imphes in turn that an ef- 
fective interaction should reflect such couplings with the 
continuum, i.e., a consistent many-body scheme should 
include bound states, resonances and virtual states as 
well. 

The present approach may allow for such a scheme. 
The simplest extension is to consider in-medium scatter- 
ing of two nucleons in e.g., infinite nuclear matter or neu- 
tron star matter, as done by Dickhoff and et al. [s^l or 
Schulze et al. [Sj . This is easily accomplished by insert- 
ing an exclusion operator Q in Eq. for the Green's 
operator. This exclusion operator can be constructed so 
as to prevent scattering into occupied states. Typically, 
one can then generate selected hole-hole, particle-hole 
and particle-particle intermediate state correlations that 
reflect a given nuclear medium. Combined with a self- 
consistent determination of the single-particle energies it 
is then possible to derive various classes of diagrams at 
the two-body level. Such work is in progress, based on 
the present method and the exact definition of the exclu- 
sion operator Q by Suzuki et al. j32| . 
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